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O . ABSTRACT 

We study the conversion of fast magneto-acoustic waves to Alfven waves by 
means of 2.5D numerical simulations in a sunspot-like magnetic configuration. A 
' fast, essentially acoustic, wave of a given frequency and wave number is generated 

^ ■ below the surface and propagates upward though the Alfven/acoustic equiparti- 

• tion layer where it splits into upgoing slow (acoustic) and fast (magnetic) waves. 

^ ' The fast wave quickly reflects off the steep Alfven speed gradient, but around and 

^ ■ above this reflection height it partially converts to Alfven waves, depending on 

the local relative inclinations of the background magnetic fleld and the wavevec- 
tor. To measure the efficiency of this conversion to Alfven waves we calculate 
acoustic and magnetic energy fluxes. The particular amplitude and phase rela- 
tions between the magnetic fleld and velocity oscillations help us to demonstrate 
^ ' that the waves produced are indeed Alfven waves. We flnd that the conversion 

^ ■ to Alfven waves is particularly important for strongly inclined flelds like those 

existing in sunspot penumbrae. Equally important is the magnetic fleld orien- 
tation with respect to the vertical plane of wave propagation, which we refer to 
as "fleld azimuth". For fleld azimuth less than 90° the generated Alfven waves 
continue upwards, but above 90° downgoing Alfven waves are preferentially pro- 
^ ■ duced. This yields negative Alfven energy flux for azimuths between 90° and 

■ 180°. Alfven energy fluxes may be comparable to or exceed acoustic fluxes, de- 

pending upon geometry, though computational exigencies limit their magnitude 
in our simulations. 
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C/3 



^Institute dc Astrofi'sica dc Canarias, 38205, C/ Via Lactea, s/n, La Laguna, Tcncrifc, Spain 

^Dcpartamento de Astrofi'sica, Universidad dc La Laguna, 38205, La Laguna, Tencrifc, Spain 

•^School of Mathematical Sciences and Monash Centre for Astrophysics, Monash University , Clayton, 
Victoria, 3800 , Australia 



-2- 



Introduction 



Observations using the Solar Optical Telescope (SOT) aboard Hinode (De Pontieu et al. 



unambiffu- 



20071 ) and the Coronal Multi-Channel Polarimeter (CoMP) fiTomczyk et al. 
ously reveal ubiquitous Alfvenic oscillati ons in the solar corona, with implications both for 
the Sun's atmosphere and the solar wind. iTomczyk et al.l find that coronal Alfvenic power is 
broadly spread in frequency, but with a distinct peak around 3-4 mHz, characteristic of the 
Sun's internal p-mode wavefie ld. Whether these trarisverse oscillations are strictly Alfven 
waves or instead kink waves (jVan Doorsselaere et al.l l2008l ) depends on the magnetic and 
density structuring of the atmosphere, and may vary with height as cross-field structuring 
becomes more or less important. We might expect though that sunspot atmospheres, with 
their presumed large-scale magnetic structures, present the most ideal site for more-or-less 
pure Alfven waves to propagate. 

The generation of Alfven waves at the photosphere and their prop agation through the 

vario us layers of t he solar atmosphe r e has been extensively modelled (e.g. lCranmer fc van Ballegooijen 
20051 ). However, ICally fc Hansen! ( 1201 ll ) recently suggested that Alfven waves must also 
be produced by mode conversion from fast magneto-acoustic waves beyond their refiection 
height in the low chromosphere. This coupling only occurs if wave propagation is not co- 
planar with gravity and magnetic field, and so the problem is necessarily three dimensional 
(3D). 

Conversion from fast-mode high-/3 magneto-acoustic waves (manifesting as p modes in 
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1 for a review. In a two-dimensional situation, the 



transformation from fast to slow magnetoacoustic modes is demonstrated to be particularly 
strong for a narrow range of magnetic field inclinations around 20-30 degrees to the vertical. 
For this reason, and because of the reduction in acoustic cutoff frequency afforded by strong 
inclined magnetic fields, magn etic field concentrat ions on the solar surface may truly be 
called magnetoacoustic portals ( iJefferies et al.ll2006l ). coupling the Sun's interior oscillations 
to those of its atmosphere. 

The remainder of the wave-energy flux though, which is near-total away from these 
preferred inclinations, enters the low-/? atmosphere as fast, predominantly magnetic waves. 
Due to the steep Alfven speed gradient with height, these fast waves soon reflect back 
downward at the height 2rcf at which uo ~ VAkh-, where u is the frequency, va the Alfven 
speed, and kh the horizontal wavenumber. 
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Fast-to- Alfven conversion occurs around and above this fast wave refl ection height Zref 
in 3p (see Fig. [1]), locahzed more closely to z^^f as frequency increases fICally fc Hansen 
201ll ). However, at the 3-5 mHz frequencies characteristic of p-modes, the process is typ- 
ically sprea d over much of the chro mosphere. Studies of fast-to-Alfven conversion were 
initiated by ICally fc Goossenj ( l2008l ). who found that it is most efficient for preferred field 
inclinations from vertical between 30 and 40 degrees, and azimuth angles (the angle between 
the vertical magnetic and wave propagation planes) between 60 and 80 degrees, and that 
Alfvenic fiuxes transmitted to the upper atmosphere can exceed acoustic fiuxes in some cases. 
Newington fc Callyl (120101 ) studied the conversion properties of low-frequency (~ 1-2 niHz) 
gravity waves, showing that even larger magnetic field inclinations can support gravity wave 
to Alfven conversion and resulting Alfvenic wave propagation to the upper atmosphere. In 
nonuniform magnetic field, the relevant angles for mode conversion, either fast-to-slow or 
fast-to-Alfven, are those pertaining in the respective conversion regions. Fig. [T] summarizes 
the overall picture of conversion between the fast, slow and Alfven waves. 

Motivated by these recent studies, we attack the problem by means of 2.5D numerical 
simulations. The purpose of our analysis is to calculate the efficiency of the conversion from 
fast-mode high-/? magneto-acoustic waves to Alfven and slow waves in the upper atmosphere 
in a spreading sunspot-like magnetic field configuration. Our i nitial results on the c onver sion 
to Alfven waves in simple field configurations were reported in lKhomenko fc Callyl (120 111 ). In 
the present work we extend our simulations to a more realistic case of a sunspot atmospheric 
models spanning the shallow subphotosphere (—5 Mm) to the high chromosphere (+1.9 Mm). 

Above this layer, the transition region acts as a partial Alfven wave refiector, affecting 
their distr ibution in the corona and solar wind and producing Alfven turbulence via nonlinear 
coupling (ICranmer fc van Ballegooijenl l2005l l. but this is beyond the scope of our present 
modeling. 



2. Simulation setup 



We numerically solve th e non-linear equations o f ideal MHD usin g our code Mancha 
(IKhomenko fc ColladosI l2006l : iKhomenko et al.l l2008l : iFelipe et al.l |2010|) • The code solves 
non-linear equations for perturbations, the equilibrium state being explicitly rem oved from 



the e quations. In this study, for simplicity, and coherence with previous studies (see lCally fc Hansen 



201ll ). we use a 2.5D approximation. This approximation means that we allow all vectors 
in three spatial directions, but the derivatives are taken only in two (one vertical and one 
horizontal) directions, so that our perturbations are only allowed to propagate in the X — Z 
plane. In addition, the initial perturbation is kept small to approximate the linear regime. 
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Fig. 1. — Schematic diagram illustrating the various mode conversions and reflections as 
a seismic ray (labeled "Fast Wave (acoustic)") enters the solar atmosphere in a region of 
strong inclined magnetic field. Field lines (pale blue) are oriented out of the plane and are 
shown here in projection. Their orientation is given by the inclination angle 9 from vertical 
(0 < ^ < 90°) and azimuth angle measured clockwise from the wave propagation plane. 
By symmetry we need only consider < </> < 180°, with < 90° in the left panel and 
> 90° in the right panel. First, at the Alfven- acoustic equipartition level va = cs the 
ray splits into an essentially acoustic field-guided slow wave (depicted in red) and a fast 
magnetically dominated wave (black). The slow wave may or may not reflect depending on 
whether u < UcCosO. The fast wave goes on to reflect higher in the atmosphere due to the 
rapidly increasing Alfven speed with height. On its way downward it again mode converts 
at the equipartition level. In the scenario depicted on the left (0 < 90°), the upward slow 
wave is much stronger than the downward one because the fast ray is more closely aligned 
with the magnetic field (small attack angle) on the upstroke than on the downstroke. This 
situation is reversed if the magnetic field were inclined in the opposite direction (equivalent 
to > 90°, right panel). In a nebulous region around and above the fast wave reflection point 
(which may extend far higher than the fuzzy blob used to represent it here), fast-to- Alfven 
conversion occurs, in the case on the left predominantly to an upgoing Alfven wave. For the 
> 90° (right), the downgoing Alfven wave is favored. The fast-to- Alfven conversion may 
only occur where the wave vector and the magne tic fiel d lines are not in the same vertical 
plane. This diagram extends the description in ICallyi (120071 ) by including fast-to- Alfven 
conversion. 
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As a background magnet o-static model atmosphere we use one sunspot-like model from 
Khomenko fc CoUadosI (120081 ) . For simplicity, this model is azimuthally symmetric, with no 
twist of the magnetic field lines. Whilst twist is often seen in sunspot magnetic fields, it is 
not a necessary feature of the processes we wish to explore here. It will be added along with 
other physical features in later studies. The computational region is magnetized in all its 
volume (distributed currents), but the strongest magnetic field is concentrated around the 
axis of the structure. The dimensions of the simulated domain are 78 Mm in the horizontal 
X direction and 7.4 Mm in the vertical Z direction. The axis of the sunspot is placed at 
the middle of the domain at X = 39 Mm. The bottom boundary of the domain is located 
at —5 Mm below the photospheric level, Z = 0. This zero level is taken to be the height 
where the optical depth at 500 nm is equal to unity in the quiet Sun atmosphere 39 Mm 
away from the sunspot axis. The thermo dynamic variables of the atmosphe re at 39 Mm 



from the axis are taken from Model S of IChristensen-Dalsgaard et a 



in the deep 

sub-photosphere layers and continuing according to the VAL-C model (IVernazza et al.lll98ll ) 



in the photospheric and chromospher ic layers . The sunspot axis in the atmospheric layers is 
given by the semi-empirical model of lAvrettl (jl98l[ ). The magnetic field at the axis is about 
900 G at Z = Mm. This is quite moderate for a sunspot, but is adopted for numerical 
reasons. The Alfven speed at the top of our computational domain becomes quite excessive 
if the field strength is too high, thereby necessitating impractically small time steps. Taking 
larger field strength would mostly result in a d ifferent scaling of the p roblem, without new 



physical phenomena being introduced (see, e.g. iKhomenko et al.ll200 



3_pr( 

3). 



The spatial resolution of our simulations is 150 km in the horizontal direction and 50 
km in the vertical direction (520 by 148 grid points). The upper 5 00 km (10 grid points) 



of the domain are occupied by the absorbing PML boundary layer (IKhomenko et al.l 12008 



Felipe et al.ll2010l ). so the effective physical boundary of the domain is located at 1.9 Mm 
above the photosphere. The PML, "Perfectly Matched Layer", is a numerical device that 
acts as an excellent absorber of waves, and thereby removes unphysical reflections from 
the edges of the computational domain. It is well-tested in the current code, and performs 
ad mirably. No sign ificant reflections are detected. The concept of PML was firstly introduced 
by iBerengerl (119941 ) for electromagnetic waves, but was quickly ext ended to other wave types. 
It is n ow u sed in many codes mo deling wave propagation, e.g. iParchevsky fc Kosovichev 
(120071 ) and iHanasoge et al.l (12010! ) . The calculations throughout are adiabatic and we use 
the ideal equation of state to close the system. 

As our modelling is 2.5D, we have made cuts through the sunspot model at several Y 
positions. Here we will describe three simulations situated in vertical planes aXY = 7.5, 11.25 
and 15 Mm from the sunspot axis. A cut through the axis is uninteresting, as it is strictly 
2D and can produce no Alfven waves. Figure [2] depicts variations of some characteristic 
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parameters of the Y = 7.5 Mm cut through the sunspot model. 

We drive waves by an imposed perturbation in a few grid points near the bottom 
boundary of the domain a.t Z = —5 Mm. The perturbation is calculated analytically as 
an acoustic-gravity wave of a given frequency and wavenumber, neglecting the magnetic 
field (dynamically unimportant at this depth), and neglecting the temperature gradient. We 
in troduce self-consistent pert urbations of the velocity vector, pressure, and density according 
to 



Mihalas fc Mihalad f ll984[ ): 



(5K = Voexp {^-^ + k^iZ^ sin(wt - k^rZ - k^x) 



or ( Z \ 

— = Vo\P\ exp I — + k^iZ ) sin{ut - k^rZ - k^x + (j)p) 
Fo \2H ) 

op / z \ 

— = Vo\R\ exp — + k^iZ sm{ut - k,rZ - k^x + (pn) 
Po V2±/ / 

SVx = Vo\U\ exp (^-^ + k^iZ^ sin(wt - k^rZ - k^x + (pu) , 
where the amplitudes and relative phase shifts between the perturbations are given by 



\R\ - " Jk^ +(k + t^'M _ 1^ ^3) 
|f/| = ^|P| (4) 
<t>P = ct>u = arctan + (5) 

Given the wave frequency w, and the horizontal wavenumber k^, the vertical wavenum- 
ber kz is found from the dispersion relation for acoustic-gravity waves in an isothermal 
atmosphere: 

K = hr + ikz, = ^{uj^-ujl)/cl-kl{uj^-ujl)/uj^ (7) 



where ujc = ^g/2cs is the acoustic cut-off frequency and ujg = 2ijjc\/^ — 1/7- 

In all simulations described in this work we set the perturbation frequency v = uj/2ti = 
5 mHz and horizontal wave number k^ = 1.37 Mm"^. According to our sunspot model 
(see Fig. |2]), the driving frequency is just slightly below the maximum cut-off frequency 
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reached at the temperature minimum. Over most of the spot these values result in the fast 
wave reflection level Zj-ei being higher than the Alfven/acoustic equipartition surface z^^, i.e., 
u/kx > cs (automatically satisfied for a propagating acoustic wave cu^ = c|(/i;^ + A;^)). This is 
important, since it allows upcoming acoustic (fast) waves in < cs to convert to magnetic 
(fast) waves in va > cs before reflecting at Zj-ef. 

To separate the Alfven mode from the fast and slow magneto-acoustic modes in the 
magnetically dominated atmosphere {va ^ cs) we use velocity projections onto three char- 
acteristic directions: 

eiong = [cos sin ^, sin sin ^, cos^]; 

Cperp = [— COS slu^ 9 slu 0, 1 — slu^ 9 siu^ 0, 

(8) 

— COS 9 sin 9 sin 0] ; 

etrans = [- COS 6^, 0, COS slu 6'] . 

These projections were shown to be rather efficient in separating the pe rturbations corre- 



spon ding to all three modes both for idealiz ed magnetic field co nfigurations (IKhomenko fc Gaily 



201ll ). and also for more complex ones (IFelipe et al.l |2010| ). The first projection (ei, 



long; 



selects the slow magneto-acoustic wave, propagating parallel to the field; the second one 
(cp firp) selects the Alf v en wa ve, according to the asymptotic polarization direction derived 



by ICally fc GoossensI (120081 ): the remaining orthogonal direction (ctrans) selects the fast 
magneto-acoustic wave. In the rest of the paper we will address these projections as "acous- 
tic", "Alfven" and "fast". 

To measure the efficiency of conversion to Alfven waves near and above the cs = va 
equipartition layer, we calculate the time-averaged acoustic and magnetic energy fiuxes as far 
as possible from the conversion layer in the magnetically dominated atmosphere {va ^ c^): 

fm« = {*Bx(il7xB„))/fi„. 

The positive sign of the fiuxes means energy propagating upwards. 

We also calculate a measure of the time-averaged energy contained in all three wave 
modes according to 

^long = PoCsiSV^Lg) 

Epcrp = poVA{SV^orp) (10) 



pcrp 

^trans = P0^A(5K?ans) ' 



where in each case we use the corresponding velocity projections into characteristic directions 
(Eq. [8]). For pure acoustic and Alfven waves, where there is strict equipartition between 
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Fig. 2. — Properties of the magneto-static sunspot model in the Y = 7.5 Mm shce away from 
the axis. The background image is the acoustic cut-off frequency unadjusted for magnetic 
effects. The sohd black line is the level where acoustic and Alfven speeds are equal {cs = va)- 
The dashed line is the fast mode reflection level calculated as the level where the wave 
frequency u and wave horizontal wave number kx are related by w = VAkx- The magnetic 
field lines are shown as colored lines, different colors meaning different azimuth of the field 
from green (0 ^ 0°) to red (0 ^ 90°) and yellow (0 ^ 180°); the is counted clockwise from 
the vertical wave propagation plane XZ. 

Vertical and horizontal dimensions are not to scale. 
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Fig. 3. — Variations with time (horizontal axis) and height (vertical axis) of the three 
orthogonal components of the velocity in the simulation at Y = 7.5 Mm from the sunspot 
axis. The color coding is the same in the three panels. The velocities are scaled with a 
factor of -x/pqCs (left panel) and ^/poVA (middle and right panels). Horizontal solid line is the 
level where cs = va', horizontal dashed line is the fast mode reflection level. Velocities are 
taken at one horizontal X location where the magnetic field inclination 6 = 31° and azimuth 
= 36°. While the fast mode (visible in Vtrans above the cs = va level) is reflected, the 
Alfven mode (visible in Vperp above the cs = va level) continues above. Note the difference 
in the propagation speeds of the acoustic and Alfven wave above cs = va- 





Fig. 4. — Snapshots of the three orthogonal components of the velocity taken at about 19 min 
after the start of the simulation at y = 7.5 Mm from the sunspot axis. Upper panel: V\ong', 
middle panel: Vtrans! bottom panel: Vperp. The blue-orange colors mean negative-positive 
velocity directions; the range of the color coding is the same in the three panels. The 
velocities are scaled with a factor of ^/poCs (upper panel) and ^PqVa (two lower panels). 
Horizontal solid line is the level where cg = va] horizontal dashed line is the fast mode 
reflection level. Magnetic field lines are inclined black lines. The axes are not to scale. Note 
the presence of the Alfven mode in V^perp above cs = va, where the fast mode (Krans) is 
already reflected. Since the velocity polarizations 'long', 'trans', and 'perp' are based on 
va ^ Cs asymptotics, care must be taken not to over-interpret these figures well below the 
equipartition level. The apparent "discontinuity" seen in Vp^j-p near Z = —1 Mm (seen also 
in Fig. [3]) is a node. 
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kinetic and compressional or kinetic and magnetic energies respectively, these would indeed 
be the true energies. Be that as it may, these forms are convenient for purposes of exposition 
and shall be used here. 



3. Velocity projections 

Figure [3] shows an example of the projected velocity components in our calculations as a 
function of height and time. The velocities are taken ai X = 10.5 Mm, where the distinction 
between the three wave modes in the magnetically dominated atmosphere is clearly visible. 
In this representation the larger inclination of the ridges means lower propagation speeds and 
vice versa. Note, that projecting the velocities according to Eq. (j8]) allows us to separate the 
wave modes only in the magnetically dominated atmosphere, above the horizontal solid line 
in Fig. [3l The figure shows how the incident fast mode wave propagates to the equipartition 
layer cs = va gradually changing its speed. After reaching the equipartition layer at about 
5 min into the simulation, it splits into several components. The essentially magnetic low-/3 
fast mode is produced above Z = 0.2 Mm (middle panel). This mode is reflected back down 
a few minutes after it has been produced. The reflection height, calculated as the height 
where the wave frequency u and wave horizontal wave number are related by a; = v^kr^ 
(ignoring the sound speed contribution to the fast wave speed), is well reproduced in the 
simulations, as the velocity variations associated with the fast mode decay rapidly above its 
reflection layer. 

Unlike the purely 2D case, there is an Alfven mode produced above the cs = va level 
(right panel). This mode has a clearly distinct propagation speed compared to the acoustic 
mode (left panel). The ridges of the Alfven mode are nearly vertical since the Alfven speed 
at these heights is large (above 100 kms~^). 

The low-/3 essentially acoustic slow mode escapes to the upper atmosphere tunnelling 
through the cut-off layer due to the field inclination of 6* = 30° that effectively reduces the 
cut-off frequency by the factor cos 6. 

The velocity amphtudes in Fig. [3] are scaled with a factor of ^^/PoVpl^, where w ph = cs for 
the V^ong component and Vp^ = va for the other two components. Although va ^ cs above 
the solid line in Fig. [3l the scaled amplitudes of the Alfven mode are still smaller than of 
the slow acoustic mode at the selected X location. 

Figure [3] also shows that the simulation enters into a stationary stage after about 10-15 
min. A snapshot of the wave field developed in the stationary stage of the simulations is 
given in Figure ID The movie of this simulation is available electronically. 
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The behaviour of the three wave modes changes across the sunspot radial direction, 
as expected due to the change of the atmospheric properties. The magnetic field is more 
inclined at the periphery of the sunspot (much more so than is suggested by the figures, 
since the vertical scale has been stretched relative to the horizontal scale) and the azimuth 
varies from = 0° to 180° from the right to the left hand side of the sunspot. The cut-off 
frequency reaches its maximum of 5.02 mHz at X = 0, Z = Mm, decreasing steeply above 
this due to the rise of temperature in the chromosphere (Fig. |2]). 

This steep temperature increase produces a partial reflection of the slow acoustic mode in 
the magnetically dominated atmosphere. This reflection is visible as an interference pattern 
formed above the solid line in the l^iong projection in the left part of the sunspot (upper 
panel). This additional reflection also produces stronger slow magnetic modes below the 
Cs = va line in the left part of the sunspot due to the secondary mode transformation from 
the downgoing waves. 

The fast magnetic mode reflection above Cg = va is complete in the central part of 
the sunspot, though the evanescent tail (the exponentially decaying part of the fast wave 
above the classical reflection point) still does not fully fit in the computational box even 
here. At the periphery, the height of the simulation domain is insufficient for the fast mode 
to complete its reflection and it is partly absorbed by the PML boundary condition. This 
produces artificially large fast mode amplitudes at the periphery of the sunspot. 

The Alfven mode component generated after the transformation is stronger at the hor- 
izontal locations where the fast wave is reflecte d at lower heigh t s (X = 0-20 Mm). The 



reason for a this is apparent from Fig. 10b of ICally fc Hansen! (120111 ): the fast-to-Alfven 
conversion region is distributed throughout the evanescent tail of the fast wave beyond its 
reflection point. At the frequencies and wavenumbers typical of our simulations here, that 
amounts to 20 scale heights or more! Since our computational box is truncated only a few 
scale heights above z^ei, we can sample only a small part of the conversion process. Where 
Zref is lower (in the center of Fig. 3c) the Alfven conversion is of course enhanced relative to 
where it is higher simply because our box includes more of the conversion region there. Our 
experiments at 25 mHz, with kr^ also increased by a factor of five, confirm this explanation, as 
for those waves the evanescent tail is much shorter and the conversion region therefore more 
compact. As expected, the higher frequency simulations yield substantially higher Alfven 
fluxes (see Section [6]). 
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4. Energy of three wave modes 



We proceed by measuring the energies contained in each of the projected velocity com- 
ponents at the upper part of our simulation domain, above the Alfven-acoustic equipartition 
layer cs = v^- As the properties of the atmosphere change in the horizontal direction across 
the sunspot, the height of the cs = va layer changes as well (see Fig. [2]). To be consistent, 
we take time averaged energies at heights from 400 km above the cs = va layer up to the 
upper boundary of our simulation box in the stationary stage of the simulations. 

Figure [5] illustrates the energies as a function of inclination and azimuth of the sunspot 
magnetic field lines at the corresponding horizontal locations. The orientation of the mag- 
netic field is taken at heights where cs = va in each of the simulations at three Y locations. 
The three strips in each panel of Fig. [5] are for the t hree simulation runs. The format of the 
figure allows comparison with the previous studies, ICally fc GoossensI ( 120081 ) Figure 2 and 
Khomenko fc Callvl feoilh Figure 4. 



This figure demonstrates that the maximum energy of slow acoustic waves is transmit- 
ted for inclinations around 30 degrees (left panel). The dependence on the inclination is 
stronger than the dependence on the azimuth. This result is consistent with previous 2D 



theoretical models of fast-to-slow mode trans 



netic field (Crouch fc Callv 



20031: 



Callv 



20061: 



brmation in the homog eneous inclined mag- 



Schunker fc Callyll2006l ). It is also consistent 



with the 3D analysis by Cally fc GoossensI ( 20081 ). The explanation of this effect is offered by 



the ray perspective. The fast-mode high-/? waves launched from the sub-photospheric layers 
(with an angle about 90 degrees, i.e. their lower turning point) reach the Alfven- acoustic 
equipartition layer with an angle close to 20 — 30 degrees (ISchunker fc Callyll2006l ). Thus, 
for the moderately inclined magnetic field the attack angle between the wavevector and the 
magnetic field lines is small and the transformation is efficient. Note also that the frequency 
of waves is just at the cut-off frequency of the sunspot atmosphere, so due to tunneling 
effects the slow mode acoustic energy is also present over the wide range of inclinations. 

The maximum energy of the Alfven wave is present at inclinations about 35° and azimuth 
between 70° and 100° (right panel of Fig. [5]). The maximum power increases from the outer 
to the inner strip indicating a tendency towards larger transmission at larger inclinations. 
Thus, the behaviour of waves in sunsp ot models is similar to the one found previously 



in the models with h omogeneous field (ICally fc GoossensI l2008l : iKhomenko fc Callyl 12011 



Cally fc Hansenll201ll ). The transmission of Alfven waves is more efficient for inclined fields, 
which is important in the sunspot penumbra. 

Some transmitted energy is also present in the -Etrans component (middle panel of Fig. [5]) 
near the sunspot periphery, though we believe this to be an artefact of our computational 
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Fig. 5. — Wave energies calculated according to Eq. [TO]at the top of the atmosphere (averages 
at heights from 0.4 Mm above the Cs = layer up to 2 = 1.9 Mm) for the three projected 
velocity components in turn, selecting respectively acoustic, fast, and Alfven waves. The 
units of the color coding are 10^ erg cm~^ s~^. The three strips in each panel are the results 
of the three simulation runs at Y = 7.5 (lower), 11.25 (middle) and 15 Mm (upper) from the 
sunspot axis. 




50 100 150 50 100 150 50 100 150 

(6 [degr] [degr] [degr] 



Fig. 6. — Logio of the amplitude ratio 5B / ^JJlqPq to 5V , both projected into the characteristic 
directions according to Eq. [HI Left panel: slow acoustic mode (eiong); middle panel: fast 
magnetic mode (ctran); right panel: Alfven mode (cperp). The three strips in each panel 
are the results of the three simulation runs at y = 7.5 (lower), 11.25 (middle) and 15 Mm 
(upper) from the sunspot axis. 
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Fig. 7. — Upper panel: phase shift between the variations of SVperp and ^Sperp projections. 
The format is the same as Fig. [61 Lower paneh same quantity, but as a function of height 
and horizontal distance in the simulation at y = 7.5 Mm from the axis. The phase shift is 
only shown above the Alfven-acoustic equipartition layer, where the projections allow the 
separation between the modes and where the Alfven mode exists after the transformation. 
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box not being tall enough to allow complete fast wave reflection there. 



5. Alfven mode polarization relations 

To check if projecting the velocities in the directions given by Eq. gives an efficient 
way to separate the Alfven mode from the fast and slow magneto-acoustic modes, we make 
use of the polarization relations for the Alfven wave. In a classical Al fven wave variations of 



the magnetic field and velocity should be related as (see |Priestlll984l ) 



SV = T^. (11) 

where the upper sign is for Alfven waves propagating parallel to the magnetic field and 
the lower sign is for those propagating in the opposite direction. The kinetic and magnetic 
energy for an Alfven wave should be in equipartition, so the ratio R = 5B / ^/]IqPq/5V should 
be equal to one. 

To confirm the Alfven nature of the transformed waves, we checked the amplitude 
and phase relations for all three modes reaching the upper atmosphere. Figure |6] presents 
calculations of the amplitude ratio R for the three modes. In each of the cases 5B and 5V 
pairs are the projections in the corresponding characteristic direction for each mode (Eq. |H]). 
The ratio is taken at the same heights as the energies from Fig.|5]and is averaged in time in the 
stationary stage of the simulations. The ratio turns out to be different for the three modes. 
In the case of the slow acoustic mode (left panel), the magnetic field variations associated 
with the velocity variations are very small (ratio R about 10~^). In contrast, in the case of 
the fast magnetic mode, magnetic field amplitudes are relatively strong providing the ratio 
about 10^ (middle panel). For the Alfven mode the ratio stays around 10° (right panel) 
meaning that the velocity and magnetic field oscillations are in equipartition as expected for 
an Alfven wave. 

Figure [7] shows the phase shift between the 5i?perp and 5V^erp corresponding to the Alfven 
mode. The situation here is very different from what one might naively expect. It is also 
different to the case of the simple atrnosph ere and homogeneous magnetic field considered by 



us previously in lKhomenko fc Callyi (1201 if ). In the right part of the sunspot, for the azimuth 
between 0° and 90° the Alfven waves mostly behave as they should. The phase shift is about 
180°, meaning an upward wave propagation. Nevertheless, in the left part of the sunspot 
(0 between 90° and 180°) the value of the phase shift indicates a downward propagation. 
Taking into account that the amplitude ratio (Fig. [6]) gives clear evidence for the Alfvenic 
nature of the waves, we conclude that for %p between 90° and 180° downward-propagating 
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Vertical component of the energy fluxes calculated after Eq. [91 Upper panel: 



pcrp) 



acoustic flux; Middle panel: magnetic flux; Bottom panel: magnetic flux, but using 5B 
^V^perp in Eq- El i.e. magnetic flux due to Alven waves. In the two bottom panel only values 
above the Alfven-acoustic equipartition layer are shown where they are meaningful. Red 
colors mean upward flux, blue colors mean downward flux. 
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Alfven waves are generated by the mode transformation. 



This accords perfectly with the uniform field modeling results of ICally fc HansenI fl201ll ) 
(see their Figures 4 and 5) indicating that between 0° and 90° the upward propagating fast 
waves couple most efficiently to upward Alfven waves. Since an alignment is needed between 
the direction of the wave propagation and magnetic field for efficient coupling, where the 
azimuth is between 90° and 180° the strongest coupling happens between the downward 
propagating fast waves (after their reflection in the magnetically dominated atmosphere) 
and the Alfven waves. In this case downward propagating Alfven waves are produced. We 
believe the same happens in the simulations, producing downward propagating Alfven waves 
in the left part of the sunspot for ip above 90°. 

All in all these calculations conflrm the Alfven nature of waves selected by the Cperp 
projection in our simulations, including returning the predicted result regarding the coupling 
between refracted fast waves and downward propagating Alfven waves in regions where the 
magnetic alignment is favourable. 



6. Energy fluxes 

We proceed by calculating the acoustic and magnetic fluxes according to Eq. 1^. Fig- 
ure M gives the time averages of the vertical magnetic and acoustic fluxes over the upper 
part of the simulation domain. In the case of the magnetic fluxes we only show them above 
the Alfven- acoustic equipartition layer where they are produced and where the separation 
between the modes according to Eq. (jS]) is meaningful. 

Above the equipartition layer, the upward acoustic flux is mostly present in the right part 
of the sunspot for intermediate fleld inclinations (upper panel). This result is in agreement 
with the energy calculations from Fig. |5] where the maximum of the transmitted acoustic 
energy is for inclinations about 30°. In the left part of the sunspot the acoustic flux is about 
zero or slightly negative. A large negative flux is also present at the central part around the 
axis at the very top of the simulation domain. Thus, slow acoustic waves are reflected in the 
left part of the sunspot. This appears to be in agreement with the visual impression from 
Fig. H] (upper panel) where the interference pattern of the slow acoustic modes is observed 
to the left of the axis. We believe that this partial reflection is due to steep temperature 
increase around the axis in the upper chromospheric part of the simulation domain. We have 
performed additional simulations placing the upper boundary of the domain at lower heights, 
below the chromospheric temperature increase. In this case we did not observe downgoing 
slow acoustic waves. This conflrms that the reflection is a physical effect and not a numerical 
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Fig. 9. — Similar to Fig. [HI but for a 25 mHz wave with = 6.85 Mm~^. The higher 
wavenumber shortens the evanescent tail of the fast wave and therefore the Alfven conversion 
region is more compact and fits within the computational box. A narrower box and three 
times finer horizontal resolution has been used here to numerically resolve short-period waves. 
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artefact from the upper boundary condition of the simulations. 

The magnetic flux above the Alfven-acoustic equipartition layer (middle panel of Fig. [S]) 
is mostly positive. Note that the absolute value of the magnetic flux is at all locations lower 
than the acoustic flux. The exception is the peripheral regions of the sunspot model where 
relatively large magnetic flux is obtained. However, this artificially large flux is due to 
insufficient height of the simulation box above the equipartition layer making it impossible 
to complete the reflection for the fast magnetic modes. 

The magnetic flux in the middle panel of Fig. |8] is due to a mixture of the fast and 
Alfven modes. To separate the flux of the Alfven mode we calculated the magnetic flux from 
Eq. (|9]) using the projections 5-Bperp and ^l^perp- The result is given in the bottom panel of 
Fig. IHl In agreement with previous considerations we observe that there is no signiflcant 
flux in the peripheral part of the sunspot, conflrming its origin due to incomplete reflection 
of the fast magnetic mode. As was already clear from the calculation of the phase relations 
(Fig. [7]) the Alfven flux in the right part of the sunspot {ip between 0° and 90°) is positive 
and the flux in the left part of the sunspot {ip between 90° and 180°) is negative, indicating 
downward propagating Alfven waves. These downward propagating waves are generated 
due to the coupling between the downgoing fast magnetic mode (after its reflection) and the 
Alfven mode. 

The Alfven flux only represents a minor fraction of the acoustic flux in Fig. [Si However, 
as explained at the end of Section [31 we believe this to be an artefact of the truncation of the 
broad Alfven conversion region. To overcome this limitation, we performed an experiment 
with flve times larger frequency {u = 25 mHz) and horizontal wavenumber {kx = 6.85 Mm~^) 
allowing us to signiflcantly compact the conversion region. This high frequency experiment 
is illustrated in Fig. [HI where indeed the upward Alfvenic flux in the right half of the region 
now greatly exceeds the acoustic flux, and there is signiflcant but very compact downward 
Alfven flux in the left half. Thus, the conclusion is that Alfven flux can actually exceed 
acoustic flux. Besides, Alfven waves are more able to propagate to great heights than are 
acoustic or fast waves, which are limited respectively by shock formation and by reflection. 
This suggests that they do indeed represent a signiflcant product of wave propagation and 
conversion in sunspots. 

7. Conclusions 

In this paper we have investigated the efficiency of conversion from fast to Alfven waves 
by means of 2.5D numerical simulations in a complex sunspot-like magnetic field configu- 
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ration. It is important to realize that quantitatively simulating mode transformation nu- 
merically is a challenge, as any numerical inaccuracies are amplified in such second-order 
quantities as wave energy fiuxes. The tests presented in this paper prove the robustness of 
our numerical procedure and offer an effective way to separate the Alfven from magneto- 
acoustic modes in numerical simulations. 

In general, the conclusions from the previous models in simplified uniform magnetic 
field configurations (ICallv fc Goossenslboosi iKhomenko fc Callvllioill i lCallv fc Hansenlboill ) 
apply also to our more complex sunspot-like magnetic field configuration, though there are 
some apparent differences. We note three points in particular. 

Firstly, the maximum of the magnetic energy of Alfven waves transmitted to the upper 
atmosphere is shifted toward more inclined fields compared to the homogeneous field case. 

Secondly, the amount of magnetic energy due to Alfven waves transmitted to the upper 
atmosphere is about 10 times lower than for the acoustic waves. This differs from the 
conclusions reached by lCally fc Goossend (120081 ) who find that at some particular inclinations 
and azimuth angles the magnetic fiux is larger. This discrepancy is entirely to be expected, 
and is an artefa c t of t he limitations of our present simulations. As seen in Figure 10b of 
Gaily fc HansenI ( 120111 ). the fast-to- Alfven conversion region for frequencies comparable to 
our 5 mHz is spread more or less uniformly over some 20 scale heights from Zrei upwards. 
Glearly, we do not have the luxury of such abundant space in our simulation box, and so 
only a small fraction of the total potential Alfven fiux is produced. In theory, this could be 
remedied in any of several ways: 



1. The box could be made substantially taller. However, this would involve encompassing 
yet larger Alfven speeds, especially if a transition region is included, with unfortunate 
consequences for our numerical time-step and therefore for the practical feasibility of 
the calculation. 

2. The magnetic field could be increased in magnitude. This is an attractive course of 
action in any case, as the 900 G at Z = adopted here on the axis is very conservative, 
being a factor of 3 less than might be expected in a mature large sunspot. This would 
lower both the cs = va equipartition surface and the u = wa^x fast refiection level 
opening up more space above for mode conversion. Of course, the maximum Alfven 
speed would again increase, with implications for the numerical time step. 



3. 



Higher hor izontal wavenurn b er wa ves could be modelled, increased by say a factor of 5. 
As seen in lGally fc HansenI (120 111 ). Figure 10b, conversion of waves with k, = k^H > 1 
(where H is the density scale height) occurs in a considerably more compact region 



- 21 - 



near 2;ref. If frequency is increased by the same factor, then u/k^ and therefore z^ei is 
unchanged, and the Alfven conversion region should fit within our computational box. 
Of course, little power is expected at such high frequencies in the Sun, but nevertheless 
it is a useful experiment to prove the point. Results of simulations on a similar model 
but with kx = 6.85 Mm~^ and a frequency of 25 mHz are shown in Fig. |9l and are in 
accord with expectations, yielding far higher Alfven fluxes. 



Finally, where the reflecting fast waves meet similarly inclined magnetic field in the 
conversion region as they propagate upward, they preferentially transfer their energy to 
upgoing Alfven waves. However, if their upward path crosses field lines at large angles there 
is little transfer. But then this correspondence may be achieved on their way back down, in 
which case downward propagating Alfven waves are the beneficiaries of their largess. This 
more efficien t coupling to downward waves happens for azimuth angles above 90°, and was 
predicted by ICally &: HansenI ( l201ll ) based on their cold plasma {/3 = 0) analysis. 



In our future studies we will extend this analysis to fully 3D simulations and include 
both magnetic field twist and a chromosphere-corona transition region. Of course twist 
contributes to the field orientation in the conversion regions, and therefore to the strength of 
mode conversion, both fast-to-slow and fast-to-Alfven, but we do not anticipate significant 
novel effects beyond this. On the other hand, a transition region (TR) is expected to be 
highly reflective to Alfven waves, and therefore to greatly re duce coronal Alfven flu xes from 
chromospheric mode conversion. Results in a zero-/? model (IHansen fc Callyll201l[ ) conflrm 
this, though they indicate that substantial coronal Alfven fluxes are still possible if the fast 
wave reflection point is within a few chromospheric scale heights of the TR, or indeed if the 
fast wave fails to reflect before it reaches the TR. It will be interesting to test this in realistic 
sunspot models. 

Ultimately, our simulations will inform and be judged a gainst observations of M HD 
waves in solar magnet ic structures. The coronal observations of lDe Pontieu et al.l ( 120071 ) and 
Tomczyk et al.l (120071 ) are useful indicators of power reaching that height, though they are at 
the limits of our technology and subject to line-of-sight integration ambiguities. Wave obser- 
vations at photospheric and chromospheric levels may be more reliable and precise. In partic- 
ular, correlations between velocity and magnetic pertur bation phases are c r ucial in disentan 



gling the various MHD modes in the low atrnosphe re (INorton et al.ll200ll : iKhomenko et al. 



2003 



Settele et al.ll2002l : iFujimura &: Tsunetall2009l ). though this is a very difficult task even 



in full 3D simulations. 

Significant progress has recently been made in the identification of wave modes observed 
in sunspots and other solar magnetic features. On the one hand, simultaneous photospheric 
and chromospheric spectropolarimetric observations (such as those done with the TIP instru- 
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ment on Tenerife, Canary Islands; ICoUados et al.ll2007 ) i ntroduce the possib ility of "follow- 
ing" wave propagation with height fICenteno et al.ll2009l : iFelipe et al.ll2010b[) and obtainiii g 
phase shifts between oscillations of magnetic field and velocity ( jBellot Rubio et al. 2000 ). 



On t he other hand, bidimensional spectrometers (as e.g. IBIS instrument; ICavallini et al. 



2000 



or DOT telescope on La Palma, Canary Islands) allow us to obtain high-resolution 
2-dimensional velocity fields at different heights. Such simultaneous observations are crucial 
for our understanding of the physics of wave propagation in active regions. In sunspot um- 
brae, where the magnetic field is predominantly vertical, slow field-aligned magneto-acoustic 
waves are firmly dete cted and even reproduced in simulations including their particular ob- 
served wave pattern ( ICenteno et al.l l2006l : iFelipe et al.l l201l[ ) . In regions with an inclined 
magnetic field (sunspot penumbra and networ k canopies), a mixture of upward and down- 
ward propagating waves is often discovered (IBraun fc Lindseyl l2000l : IVecchio et al.l 120071 : 



Kontogiannis et al.ll2010[ ). At photospheric level, these up- and downgoing waves are, pos- 



sibly, fast magneto-acoustic waves undergoing a reflection p rocess. In fact, fr o m me asured 
phase shifts between velocity and magnetic field oscillations, iKhomenko et al.l (120031 ) found 
that the contribution of slow mode waves is larger in umbral regions and the contribu- 
tion of fast mode waves becomes progressively more important towar d the umbra-penumbra 
boundary. Alfven wa ves are usually detected higher up in the corona ( iDe Pontieu et al.ll2007l : 



Tomczyk et al.ll2007l ). The simulations presented here suggest that the complete decoupling 



between the fast and Alfven waves happens in the upper chromosphere, or even above. Since 
the upcoming acoustic waves on the Sun are uncorrelated, their attack angle with respect 
to the sunspot magnetic field can be arbitrary, so a mixture of up- and downgoing Alfven 
waves will be produced, mostly above regions with inclined magnetic field. Time variations 
of the chromospheric magnetic fields, together with velocity oscillations, are strongly de- 
sirable for the firm detection of Alfven wa ves. However, measuring magnetic f ields in the 
chromosphere may be challenging (see, e.g.. lManso Sainz &: Trujillo Buendl2010l ) and actual 
measurements are scar ce. As of today, very few measurements exist in active and quiet chro- 



mospheric regions (e . g. ISocas-Navarro 



2005 



Truiillo Bueno et al.ll2005l : ISocas-Navarroll2007l : 



Centeno et al.ll2010t IStepan fc Trujillo Buendl2010l ). But, as far as we are aware, no time 



variations have been reported. Further instrumental efforts, developments of chromospheric 
diagnostic techniques, as well as improved modelling, will be needed to constrain the results 
of the presented study. 
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